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■ Abstract. In this paper we consider limit theorems, symmetry of distribution, and absorption 
problems for two types of one-dimensional quantum random walks determined by 2 x 2 unitary 
matrices using our PQRS method. The one type was introduced by Gudder in 1988, and the other 
type was studied intensively by Ambainis et al. in 2001. The difference between both types of 
^sO ■ quantum random walks is also clarified. 
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1 Introduction 

The classical random walk (CRW) in one dimension is the motion of a particle located on the 
set of integers. The particle moves at each step either one unit to the left with probability 
p or one unit to the right with probability q — 1 — p. The directions of different steps are 
independent of each other. This CRW is often called the Bernoulli random walk. In the 
present paper, we consider quantum variations of the Bernoulli random walk and refer to 
such processes as quantum random walks (QRWs) here. Sometimes the QRW is also called 
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. the quantum walk. 

Very recently, considerable work has been done on discrete-time and continuous-time 
QRWs by a number of groups in connection with quantum computing. Examples include 
Aharonov et al. (2001), Ambainis et al. (2001), Bach et al. (2002), Brun, Carteret and 
Ambainis (2002a, 2002b, 2002c), Childs, Farhi and Gutmann (2002), Childs et al. (2002), 
Diir et al. (2002), Kempe (2002), Kendon and Tregenna (2002a, 2002b), Konno (2002a, 
2002b), Konno, Namiki and Soshi (2002), Konno, Namiki, Soshi and Sudbury (2003), Leroux 
(2002), Mackay et al. (2002), Moore and Russell (2001), Severini (2002a, 2002b), Travaglione 
and Milburn (2002), Yamasaki, Kobayashi and Imai (2002). For a more general setting 
including quantum cellular automata, see Meyer (1996). The present paper is concerned 
only with the discrete-time case. 

QRWs behave quite differently from CRWs. It is well known that the probability dis- 
tribution of a CRW is given by the form of a binomial distribution. The variance of the 

Key words and phrases. Quantum random walk, the Hadamard walk, limit theorems, absorption prob- 
lems. 
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CRW increases linearly with the number of time steps. By contrast, the probability distri- 
bution of a QRW has a complicated oscillatry form and the variance of the QRW increases 
quadratically with the number of time steps (see Ambainis et al. (2001), Konno (2002a), for 
examples). Such behavior of the QRW is due to the interference between the separate paths 
of the walk. 

The QRW was introduced by Gudder (see Section 7.4 in his book (1988)) to describe 
the motion of a quantum object in discrete space-time (very recently, I happened to know 
this fact from Severini (2002a)). QRWs have been re-discovered and studied intensively by 
Ambainis et al. (2001) in the context of quantum computation recently. So for simplicity we 
call the former QRW the G-type QRW, and the latter QRW the A-type QRW respectively 
in this paper. Both QRWs are essentially same, however they have a difference. One of 
the main purpose of this paper is to clarify the different behaviour between the two types 
QRWs, for examples, limit distribution, symmetry of distribution, absorption probability. 

Here we present not only our recent results on limit theorems, symmetry of distribution, 
and absorption probabilities for A-type QRWs, but also the related results for G-type QRWs, 
by using our PQRS method based on combinatorics which is different from Fourier analysis. 
This paper is organized as follows. Section 2 treats the definitions of both A-type and G-type 
QRWs and explains the PQRS method. In Section 3, we give a new type of limit theorems 
for both types of QRWs. Section 4 is devoted to absorption problems for QRWs. 

2 Definition and PQRS Method 

The time evolution of the one-dimensional QRW is given by the following unitary matrix: 

_ a b 
c d 

where a,b,c,d G C and C is the set of complex numbers. So we have 

\ a \ 2 + \b\ 2 = |c| 2 + \d\ 2 = 1, ac + bd = 0, 
c = —Ab, d = Aa 

where z is a complex conjugate of z G C and A = det U = ad — be. We should note that the 
unitarity of U gives |A| = 1. 

For the A-type QRW, each coin performs the evolution 

\L) -> U\L) = a\L) + c\R), 
\R) -> U\R) = b\L) +d\R) 

at each time step for which that coin is active, where R and L can be respectively thought 
of as the heads and tails states of the coin, or equivalently as an internal chirality state of 
the particle. The value of the coin controls the direction in which the particle moves. When 
the coin shows L, the particle moves one unit to the left, when it shows R, it moves one unit 
to the right. (The G-type QRW is also interpreted in a similar way. We will give a precise 
definition of both types of QRWs later.) In this meaning, the QRW can be considered as 
a quantum version of the CRW with an additional degree of freedom called the chirality 
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which takes values left and right. As for ways to regain the CRW from the QRW, see Brun, 
Carteret and Ambainis (2002a, 2002b, 2002c). 

The amplitude of the location of the particle is defined by a 2-vector e C 2 at each 
location at any time n. The probability the particle is at location k is given by the square 
of the modulus of the vector at k. For the j-type QRW (j = A,G), let \^j t k{n)) denote the 
amplitude at time n at location k where 



with the chirality being left (upper component) or right (lower component). Then the 
dynamics for \^ A,k{ n )) for the A-type QRW is given by the following transformation: 



(2.1) 
where 



|*A,fc(n + 1)) = P A \^ A,k+i{n)) + Q A \^A,k-i(n)) 



a b 




Qj 





c d 



It is noted that U = Pa + Qa- The unitarity of U ensures that the amplitude always defines 
a probability distribution for the location. 

On the other hand, the G-type of QRW can be determined by 



where 



\*G,k(n + 1)) = P G |* Gl fc+iW) + Qa\^G,k-i(n)) 



Pg 
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In this G-type case also, we see that U — Pq + Qg, an d the unitarity of U ensures that the 
amplitude always defines a probability distribution for the location. In our one-dimensional 
setting, Gudder considered the following simple model (see Eqs. (7.33) and (7.34) in page 279 
of Gudder (1988)) and computed \^G,k( n )) O 11 our notation) for this model (see Corollary 
7.24 in page 285 of his book) by using a Fourier analysis which is different from our PQRS 
method: 



U 



a ib 
ib a 



where < a < 1 and b = yl — a 2 . 

The simplest and well-studied example of an A-type QRW is the Hadamard walk whose 
unitary matrix U is defined by 



V2 



1 1 
1 -1 
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The dynamics of this walk corresponds to that of the symmetric CRW. In general, the 
following unitary matrices can also lead to symmetric QRWs: 



i(4>+il>) 



-V>) 



e i(4>-ip) _ e -i(0+V>) 



where r\, 0, and ip are real numbers (see pp. 175-176 in Nielsen and Chuang (2000), for 
example). In particular, we see t/o,o,o = H. 

However symmetry of the Hadamard walk depends heavily on the initial qubit state, see 
Konno, Namiki and Soshi (2002). Another generalization of the Hadamard walk is: 



H(P) 



Vp V 1 - p 
V 1 - p ~Vp 



where < p < 1. Note that p = 1/2 is the Hadamard walk, that is, H = H{l/2). 

In the present paper, the study on the dependence of some important properties and 
quantities (e.g., symmetry of distribution, limit distribution, absorption probability) on ini- 
tial qubit state is one of the essential parts, so we define the collection of initial qubit states 
as follows: 

$ = {y? = [*] GC 2 :M 2 + |/3| 2 = lj 

Let X^ n be the j-type QRW at time n starting from initial qubit state f G <£> with 
X^ = for j = A, G. In our treatment of QRWs, as well as the matrices Pj and Qj, it is 
convenient to introduce 

Ra = 



and 



Rg 
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d ' 
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, Sg = 
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We should remark that both A- type and G-type P j: Qj, Rj, Sj (j = A, G) form an orthonor- 
mal basis of the Hilbert space M 2 (C) which is the vector space of complex 2x2 matrices 
with respect to the trace inner product (A\B) = tr (A*B). Therefore we can express any 
2x2 matrix X conveniently in the form, 



(2.2) 



X = tr{P*X)Pj + ti(Q*X)Q j + ti(R*X)R j + tr(S*X)S j 



for each j = A, G. We call the analysis based on Pj, Qj, Rj, Sj (j = A, G) the PQRS method. 

The n x n unit and zero matrices are written I n and O n respectively. For instance, if 
X = I 2 , then Eq. (2.2) gives 



(2.3) 



J 2 = aPj + dQ~ + cRj + bSj 



for each j = A, G. The next table of products of Pj, Qj, Rj, Sj (j = A, G) is very useful in 
computing some quantities: 
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Qj 


Rj 


Sj 


Pj 


aPj 


bRj 


aRj 


bP 3 


Qj 


cS 3 


dQj 


cQj 


dSj 




C Pj 


dRj 


cRj 


dPj 


S 3 


aSj 


bQj 


aQj 


bSj 



where PjQj = bRj, for example. We should remark that the algebraic structure for both 
types is the same. 

In order to consider absorption problems stated in Section 4, now we describe the evo- 
lution and measurement of the A- type QRW starting from location k on {0, 1, . . . , N} with 
absorbing boundaries (for examples, see Ambainis et al. (2001), Bach et al. (2002), and 
Kempe (2002) for more detailed information). As for the G-type QRW, we can define its 
evolution and measurement in a similar fashion. 

First we consider N = oo case. In this case, an absorbing boundary is present at location 
0. The evolution mechanism is described as follows: 

Step 1. Initialize the system p e <E> at location k. 

Step 2. (a) Apply Eq. (2.1) to one step time evolution, (b) Measure the system to see 
where it is at location or not. 

Step 3. If the result of measurement revealed that the system was at location 0, then 
terminate the process, otherwise repeat step 2. 

In this setting, let (n) be the sum over possible paths for which the particle first hits 
at time n starting from k for the A-type QRW. For example, 

S£?(5) = PIQaPaQa + P\Q\ = (ab 2 c + a 2 bd)R A 

The probability that the particle first hits at time n starting from k is 

P&\n;<p) = |S<#(nM 2 

So the probability that the particle first hits starting from k for the A-type QRW is 

oo 

71=0 

Next we consider N < oo case. This case is similar to the N = oo case, except that two 
absorbing boundaries are presented at locations and N as follows: 
Step 1. Initialize the system p e $ at location k. 

Step 2. (a) Apply Eq. (2.1) to one step time evolution, (b) Measure the system to see 
where it is at location or not. (c) Measure the system to see where it is at location N or 
not. 

Step 3. If the result of either measurement revealed that the system was either at location 
or location N, then terminate the process, otherwise repeat step 2. 

Let S^f2( n ) be the sum over possible paths for which the particle first hits at time n 
starting from k before it arrives at iV for the A-type QRW. For example, 

(2.4) = P 2 a QaPaQa = ab 2 cR A 

In a similar way, we can define P A N f}(n; tp) and P A N f}(p). 
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3 Limit Theorem 

This section treats limit theorems for both A-type and G-type QRWs Xf n (j = A,G). To 
study P(Xj n = k) for n + k = even, it suffices to understand the following combinatorial 
expression. For fixed / and m with I + m = n and m — I = k, we consider 

m)= P^Q^PfQ] 12 ■ ■ ■ P l fQf n 

li,m,i>0:m\-\ \-m„=m,hA \-l n =l 

since 

P(Xl n = k) = \~ ] (Um) V \ 2 
Since Pj, Qj, Rj, Sj (j = A, G) are a basis of M 2 (C), Ej(l, m) has the following form: 



Ej(l, m) = pj(l, m)Pj + qj(l, m)Qj + rj(l, m)Rj + Sj(l, m)S. 



Next problem is to obtain explicit forms of pj(l,m),qj(l,m),rj(l,m), and Sj(l,m). In the 
above example of n = I + m = 4 case, we see that for j = A, G, 

Sj(4, 0) = a 3 Pj, Sj(3, 1) = 2abcPj + a 2 bRj + a 2 cSj, 

2, (2, 2) = bcdPj + abcQj + abdRj + bc 2 Sj, 

Zj(l, 3) = 2bcdQj + bd 2 Rj + cd 2 S j , 5,(2, 2) = d 3 Q j 

So, for example, 

Pj(S, 1) = 2abc, g J (3,l)=0, ^(3, 1) = a 2 b, Sj(3, 1) = a 2 c 

Note that Pj(l, m), qj(l, m),rj(l, m), Sj(l, m) do not depend on the type j. 

In abed = case, the argument is much easier. So from now on we focus only on abed ^ 
case. In this case, the next key lemma is obtained by a combinatorial method. 



LEMMA 1 We consider both A-type and G-type QRWs with abed ^ 0. Suppose that l,m> 
with I + m = n, then we have 
(i) for I A m{— min{/, m}) > 1, 



I Am 



Ej(l,m) = a l d rn Y, 

7=1 



|6py l- 1\ m- V 



(ii) for l(— n) > 1, m — ; 
(Hi) for I = 0, m{— n) > 1, 
where j = A, G. 



x 



cry c«7 c b 



E j (l,0) = a l - 1 P j 



Ej(0,m) = d m - Y Q 
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The proofs of parts (ii) and (iii) are trivial. The proof of part (i) is based on a consequence 
of enumerating the paths of drift I + m = n. To do so, we consider the following 4 cases: 



PjPj ' ' ' Pj QjQj ' ' ' Qj PjPj ' ' ' Pj ' ' ' QjQj 'Qj PjPj ' ' ' Pj 

W\ W2 UI3 W27 W2~,+l 



QjQj Qj PjPj ' ' ' Pj QjQj Qj PjPj ' ' ' Pj QjQj Qj 

Wl W2 W;i W2~f 



PjPj ' ' ' Pj QjQj 'Qj PjPj ' ' ' Pj ' ' ' QjQj ' ' ' Q. 



3 3 3 ^3^3 ^3 ± 3 ± 3 3 

Wl W2 W3 




QjQj ' ' ' Qj PjPj ' ' ' Pj QjQj ' ' ' Qj ' ' ' PjPj 



where w±, W2, ■ ■ ■ , u>2 7 +i > 1 an d 7 > 1- The above each type of paths of Pj and Qj 
corresponds to the each term of Pj,Qj, Rj, Sj respectively in the right hand side of the 
equation of part (i). As for the A-type case, the details of this proof appear in Konno 
(2002b). A similar proof also can be seen in Appendix A of Brun, Carteret and Ambainis 
(2002b). The proof of the G-type case is the same. 

By this lemma, the characteristic function of Xj n (j = A, G) for abed 7^ case is obtained. 
Moreover, the mth moment of X^ n can be derived from the characteristic function in the 
standard fashion. Here we give only the result of the mth moment. 



PROPOSITION 2 We consider both A-type and G-type QRWs with abed ^ 0. 
(i) When m is odd, we have 



where 



-n m Ti 



V 1 ! 1 ] k k 

+ E EE 

k=l 7=1 5=1 

<n - 2k) m+1 



\b\ 2 V +5 (k- 1\ (k- 1\ (n- k- 1\ (n- k- 1 N 



,7 — 1 1 \ 5 — lyl 7 — 1 



5-1 



x 



7<5 



= (M 2 - l&l 2 ) (l«| 2 - W) + 2(aa6/5 + aab/3), T G = \a\ 2 - \(3\ 2 , 



A A = -{n(\a\ 2 -\b\ 2 )+ 1 + 5}(\a\ 2 -\f3\ 2 ) + 

A G = (7 + 5 - n){\a\ 2 - \[3\ 2 ) - ^Qg, 
Q A = CLCib/3 + aab/3, Qq — a/3ca + a[3ca 
(ii) When m is even, we have 



'l±5_ 
. \b\ 2 



2n Q A , 



7 



M * * / |6| 2 V +< A - A A - i\ A A> - * - A 

x ("^{ (n _ fc) , + ^ 7+4)+ l7*}- 

It should be noted that when m is even, E((Xj n ) m ) is independent of initial qubit state <p 
and types j = A, G. In particular, we use the above result of m — 1 case in order to study 
symmetry of distributions for the QRWs. Moreover we have the following new type of limit 
theorem (as for the A-type QRW, see Konno (2002a)): 



THEOREM 3 We consider both_A-type and G-type QRWs with abed ^ 0. Let Q A = 
aab(3 + aab(3, and Qq = a/3ca + af3ca. If n — > oo ; then 

^ ^ ZJ (j = A,G) 
n J 

where Zj has a density 

/ ^ ,[a - fl) % ( i-^-^ { 1 -( |a| '- w ' + il) j } 

for x G ( — \a\, \a\) with 

E(Zf) = - (W - + ^ x (1 - y/l-\a\*) 

E{{zjf) = i - v/i-M 2 

and Y n =>- Y means that Y n converges in distribution to a limit Y . 



We remark that standard deviation of Zj is not independent of initial qubit state if — *[«, j3\. 
The above limit theorem suggests the following result on symmetry of distribution for both 
A-type and G-type QRWs. This is a generalization of the result give by Konno, Namiki and 
Soshi (2002) for the A-type Hadamard walk and by Konno (2002a) for the A-type QRW. 
For j — A,G, define 

$> j s = {ip g $ : P(Xf, n — k) — P(X? n = -k) for any n G Z+ and k G Z} 

*;,o = {v 9 G $ : SpTjy = for any n G Z + } 

= {v9 = '[a,/5]G$:|a| = |/3|, 6,=0} (j = A, G) 

and Z (resp. Z + ) is the set of (resp. non-negative) integers. For ip G $ s , the probability 
distribution of Xj n is symmetric for any n G Z + . Using the explicit form of E(XJ^) given 
by Proposition 2 (i) (m = 1 case), we have 
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THEOREM 4 We consider both A-type and G-type QRWs with abed ^ 0. Then for j = 
A, G we have 

$3,8 = $3,0 = 'I',. . 

From now on we give an outline of our proof of Theorem 3 (for more details, see Konno 
(2002b)). To do so, we introduce the Jacobi polynomial P^ ,M (x), where P^(x) is orthogonal 
on [—1, 1] with respect to (1 — x) u (l + x)^ with v,ji> —1. Then the following relation holds: 

(X) = T(nlt)rt+ lf Fli ~ n > n + " + ^ + 1 ^+ 1 ^ 1 - *)/2) 

where 2-^1(0, b; c; z) is the hypergeometric series and T(z) is the gamma function. In general, 
as for orthogonal polynomials, see Koekoek and Swarttouw (1996). Remark that 

5(-srK^i)(";-T 1 )-i w * ,, * r " w ' |, - i) 

t {-jff" (* : 1) ( n ~-V) = M- 2tt - v, ^w - 1) 

By using the above equations and an expression of the characteristic function for Xj n , we 

obtain the next asymptotics of characteristic function E(e t ^ x ^ n ): if n — > 00 with k/n = 
x e (-(1 - |a|)/2, (1 + \a\)/2), then 



E{e* x U n ) ~ 2 



n-l l 
2 J 

2n-4fe-2|T |4 



X 



fc=l 

,2 



2, jx + l (pfcl , n _ 2fc)2 _ + ^(po^^j cos((1 _ 2x)0 

_ (i-^)r J (p , n _ 2fc)2 _ 2/ _ ^ |2 _ w _ pM- 2fep M- 2fe |. sin((1 _ 2x)0 

where /(n) ~ means f{n)/g{n) -> 1 (n -> 00), P^ -2 * = P^l™r 2fc (2|a| 2 - 1) (i = 0, 1), 
and /a = 1,/g = -1- 

Next we prepare the following asymptotic results for the Jacobi polynomial p^+ an ^+ bn (^ x ) 
derived by Chen and Ismail (1991): if n — > 00 with k/n = x e ( — (1 — |a|)/2, (1 + |a|)/2), 
then 

9 1 „|2fc— n 

P fc iV ^ ~ P cos(An + P) 
'ixn\J— A 



where A = (1 — |a| 2 )(4a; 2 — 4x + l — |a| 2 ), A and P are some constants (which are independent 
of n), and G [0, 7r/2] is determined by cos 9 = - \a\ 2 )/4x{l-x). 
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Combining the above observations with the Riemann-Lebesgue lemma (see p. 462 in Dur- 
rett (1996)), we see that, if n — > oo, then 



12 r L 



dx 



* J1 -=^ x(l-:rW(H 2 -l)(4:r 2 -4:r+l-|a| 2 ) 



cos((l - 2x)£) - (1 - 2x) ( \a\ 2 - \(3\ 2 + i% J isin((l - 2s)0 



rial 



-(1 — x 2 )^J\a 



2 _ x 2 | \ | a 



2 



Then 0(£) is continuous at £ = 0, so the continuity theorem (see p. 99 in Durrett (1996)) im- 
plies that X%/n converges in distribution to the Z v with characteristic function <fi. Therefore 
Theorem 3 is obtained. 

Muraki (2002) introduced a notion of quasi-universal product for algebraic probability 
spaces and showed that there exist only five quasi- universal products, that is, tensor product 
(case 1), free product (case 2), Boolean product (case 3), monotone product (case 4), and 
anti-monotone product (case 5). Algebraic central limit theorems describe limit behaviors 
of rescaled sum of algebraic random variables <7i, 02, . . . with mean and variance 1, 

Q\ + cr 2 H h cr„ 

y/n 

converges weakly, in the limit n — > 00, to the Gaussian distribution e~ x2 1 2 j \/2n (case 1), the 
semi-circle distribution X[-2,2]{x)\/4: — x 2 /27r (case 2), the Bernoulli distribution (5_ 1 + 5 1 )/2 
(case 3), and arcsine distribution X(-^/2 v / 2)( ;r )/ 7r v / 2 — x 2 (cases 4 and 5), where x{A) = 1 
if x e A, = if x g A with A C R (see Muraki (2001), Hashimoto (2002), for examples). 
So our limit theorem (Theorem 3) would belong to another category which is different from 
ones given by Muraki. 

Now we compare our analytical results (given by Theorem 3) with the numerical ones 
for the A-type Haramard walk. 

We see that Theorem 3 implies that if — \/2/2 < a < b < y/2/2, then as n — ■> 00, 

P{a < X\ _/n < b — > / , — dx 

V ~ ^ ~ 1 Ja 7l(l-X 2 Wl~2x 2 

for any initial qubit state ip = t [ct,0\. For the symmetric CRW Y° starting from the origin, 
the well-known central limit theorem implies that if— oo<a<5<oo, then as n — > 00, 



J a 



e -x 2 /2 



2tt 



dx 
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This result is often called the de Moivre-Laplace theorem. When we take p> = '[l/v^i/v 7 ^] 
(symmetric case), then we have the following QRW version of the de Moivre-Laplace theorem: 
if -y/2/2 < a <b < y/2/2, then as n -> 00, 

So there is a difference between the QRW X\ n and the CRW Y° even in a symmetric case 
for ( p = t [l/y/2,i/y/5\. 

Noting that E(X^ ) = (n > 0) for any </? G $a,s, we have 

sd(X£J /n -> 7(2-^/2 = 0.54119 . . . 

where sd(X) is the standard deviation of X. This rigorous result reveals that numerical 
simulation result 3/5 = 0.6 given by Travaglione and Milburn (2002) is not so accurate. 

As in a similar way, when we take (p = *[0, e t9 ] where 9 G [0, 2ir) (asymmetric case), we 
see that if —y/2/2 < a < b < y/2/2, then as n — > 00, 

P (a <xx><^/; ^-^— * 

So we have 

E{Xl n )/n -> (2 - v^)/2 = 0.29289 . . . , sd(X\ n )/n -> \/{y/2-l)/2 = 0.45508 . . . 

When </? = *[0, 1] (^ = 0), Ambainis ei a/. (2001) gave the same result. In the paper, they took 
two approaches, that is, the Schrodinger approach and the path integral approach. However 
their result comes mainly from the Schrodinger approach by using a Fourier analysis. 

In another asymmetric case p> = t [e td ,0] where 6 G [0, 2tt), a similar argument implies 
that if -V2/2 < a <b < y/2/2, then as n — > 00, 

P{a < XlJn <b)^f 1 2 dx 

Ja 7T(1 + X)y/1 — 2x l 

Note that f A {-x; 1 [e id , 0]) = f A {x\ *[0, e id }) for any x G {-y/2/2, y/2/2). Therefore concerning 
the mth moment of the limit distribution, we have the same result as in the previous case tp = 
t [0,e te ]. So the standard deviation of the limit distribution Z\ is given by \j (a/2 — l)/2 = 
0.45508 .... Simulation result 0.4544 ± 0.0012 in Mackay et al. (2002) (their case is 6 = 0) 
is consistent with our rigorous result. 

4 Absorption Problem 

From now on we consider absorption problems for both A-type and G-type QRWs located 
on the sets {0, 1, ... , iV} or {0, 1, . . .}. Results in this section for the A-type QRW appear 
in Konno, Namiki, Soshi and Sudbury (2003). 

Before we move to a quantum case, first we describe the CRW on a finite set {0, 1, ... , N} 
with two absorption barriers at locations and N (see Doyle and Snell (1984), Grimmett 
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and Stirzaker (1992), for examples). The particle moves at each step either one unit to the 
left with probability p or one unit to the right with probability q = 1 — p until it hits one 
of the absorption barriers. The directions of different steps are independent of each other. 
The CRW starting from k G {0, 1, ... , N} at time n is denoted by Y% here. Let 

T m = min{n > : = m} 

be the time of the first visit torn G {0,1,..., N}. Using the subscript k to indicate Yq = k, 
we let 



Pf° = P k (T < Tj 



N 



be the probability that the particle hits starting from k before it arrives at N. The 
absorption problem is also known as the Gambler's ruin problem. 

Now we review some known results and conjectures on absorption problems related to 
this paper for A-type QRWs. 

In the case of U — H (the Hadamard walk), when N = oo, that is, the state space is 
{0, 1, . . .} case, Ambainis et al. (2001) proved 

(4.5) ^ ) ( t [0,l]) = P& ) ( t [l,0]) = - 

7T 

and Bach et al. (2002) showed 

for any initial qubit state (p = *[a,0\ G Furthermore, in the case of U — H(p), Bach et 
al. (2002) gave 



i im pgc[ 1 ,oi) = c03 "^- 2rt 

fe^oo 7T 

The second result was conjectured by Yamasaki, Kobayashi and Imai (2002). 

When N is finite, the following conjecture by Ambainis et al. (2001) is still open for the 
U = H case: 

Solving the above recurrence gives 

(4.6) Pf 1 ) ( t [0 i l]) = ^x ( 3 + 2 ^) Af ' 1 - 1 {N>1) 

However in contrast with Pi^(*[0, 1]) = 2/tt, Ambainis et al. (2001) proved 

lim PZ\%l]) = l/V2 
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Let Tq be the first hitting time to 0. We consider the conditional mth moment of T start- 
ing from k = 1 given the event {T < oo}, that is, E^\(T ) m \T < oo) = E^\(T ) m ;T < 

oo)/P^\T < oo) for the j-type QRW. 

From now on we begin with the classical case. In this case, to obtain P k N \ we use the 
following difference equation: 

(4.7) P W =p pW + q pW {1 < k <N-l) 
with boundary conditions: 

(4.8) P (A ° = 1, PW= 

To consider a similar equation even in the quantum case and to use the PQRS method are 
our basic strategy. 

From now on we focus on 1 < A; < iV — 1 case. So we consider only n > 1 case. Noting 
that {Pj, Qj, Rj, Sj} for each j = A, G is a basis of M 2 (C), H^(n) can be written as 

45V) = pSJ?(n)^ + Qfk\n)Qj + rg?(n)i2,- + s%\n)S 3 

From the definition of E^ k \n), it is easily shown that there exist only two types of paths, 
that is, Pj . . . Pj and Pj . . . Qj, since we consider only a hitting time to before it arrives at 
N. Therefore we see that qjj){n) = Sj N k \n) = (n > 1). 

We assume N > 3. Noting that the definition of Efl\n), we have 



4? w = 45-> - m + 45l> - w (! < k < N - !) 

The above equation is a QRW version of the difference equation, i.e., Eq. (4.7) for the CRW. 
As an example, see Eq. (2.4). Then we have 



Pj,k ( n ) = a Pj,k-i( n - 1) + crj,k-i(n ~ !) 

^•5V) = M5+i( n - x ) + rfr i5+i( n - X ) 



Note that the above equations do not depend on types of QRWs. Next we consider boundary 
conditions related to Eq. (4.8) in the classical case. When k = N, 

P${0;<p) = |S$( 0)v f = 

for any So we take S^(0) = 2 . In this case, Eq. (2.2) gives 

P$(0) = r$(0) = 

If fc = 0, then 

^(0;^) = |s(y (0)^ = 1 
for any So we choose H^(0) = J 2 . From Eq. (2.4), we have 

P$>) = a, r$>(0) = c 
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It is noted that the above boundary conditions also do not depend on types of QRWs. 
Therefore, from now on we will omit subscript j oip ( j^\n) and rj^ (n). Let 



v\ >(n) 



L A N \n) j 



Then we see that for n > 1 and 1 < k < N — 1, 



(4.9) vi N \n) 
and for 1 < k < N, 



a c 






b d 



; k+l 



n-l) 



«?>) = 



a 

c 



4 N, (o) = 








Moreover, 








(n > 1) 



So the definition of Pj^\ip) gives 
LEMMA 5 

oo 



n=i 



i^W) = Q,i(n)M 2 + C*, 2 l 



+ 2^{C j:Z {n)a(3) 



where ^(z) is the real part o/z 6 C, (p = t [oi,(3] G $ and 

C Ail (n) = \a P [ N \n) + cr{ N) (n)\ 2 , C A , 2 (n) = \h P [ N \n) + dr{ N) (n)\" 



C Afl {n) = (apf\n) + ' {n)){bpf > {n) + dr^\n)) 



Ca,i(n) = \pf\n)\\ C G , 2 (n) = \rT\n)\\ C G , 3 (n) = > '(n)^ > '(n) 



To solve P^\ip), we introduce generating functions of p^'(n) and r^ VJ (n) as follows: 



(JV), 



By Eq. (4.9), we have 



n \ z ) = J2 r l N) ( n ) zT ' 

n=l 



n=l 



V k N \z) = azpZliz) + czrf_\(z) 
rf\z)=bzpZ\(z)+dzrf + \(z) 
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Solving these, we see that both p^\z) and rj^\z) satisfy the same recurrence: 

dA%(z) - (Az + I) p[ N + \(z) + aA N \z) = 
d^)-(Az+l)rf + \(z)+a^\z)=0 

From the characteristic equations with respect to the above recurrences, we have the same 
root: if a ^ 0, then 



Az 2 + 1 =F y/AV + 2A(1 - 2\a\ 2 )z 2 + 1 
A± ~~ 2Aaz 

where A = det U = ad — be. 

From now on we consider mainly U = H (the Hadamard walk) case with N = oo. 
Remark that the definition of E^\n) gives p£ °\n) = 0(n > 2) andp^l) = 1. So we have 
Pi°°^(z) = z. Moreover noting lim^oo ffi^ (z) < oo, the following explicit form is obtained: 



jfr\z) = zX%~\ rt\z) = _1±^H X ^ 

where X± = (z 2 — l±\^z 4: + l)/\^2z. Therefore for A; = 1, we have fj 00 -* (z) = (—l+y/z 4 + l)/z. 
By using these, we obtain 

PROPOSITION 6 For each j -type QRW, we have 

(4.10) Pi^) = - + 2(l--Wa/3), 

(4.11) PS ) M = H 2 +(^-l)|/3| 2 
for any initial qubit state ip = t [a,(3] G $. 

Eqs. (4.10) and (4.11) are a generalization of Eq. (4.5) given by Ambainis et al. (2001). 
Then we get (4 — 7i)/n < P^\(p) < 1 for each j-type QRW. Moreover we obtain a result 
on the conditional mth moment of the first hitting time to starting from k — 1 given an 
event {T < oo}, that is, E$\(T ) m \T < oo) = E$\(T ) m ;T < oo)/P, ( J ) (T < oo) in a 
similar way. 

PROPOSITION 7 For each j-type QRW, we have 

E^(T \T <oo) = -J— 

n(oo). 



E%\(T ) m \T <oo) = oo (m>2) 

Next we consider a finite N case. By using boundary conditions: p[ N \z) = z and 
^n-i( z ) = 0, we see that p^\z) and rjf \z) satisfy 



(4.12) pl N \z) = g + E z ) A*" 1 + (| - £,) A*" 

(4.13) rf } (z) = C^A*"^ 1 - A^ 1 ) 
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where 

(4.14) C Z = ^=(-1) N - 2 (X^ 3 -X N - 3 ) 

x |(A^~ 2 - A^ 2 ) 2 - ^|(A^ 2 - A^ 2 )(A^ 3 - A^ 3 ) - (-1)^ 3 (A + - A_) 2 



(4.15) E z = 



N-2 \N-2\ 



2(\ 1 ;- 2 - A 



2(-l) JV - 3 (A + -A_)(Af- 3 -A^ 3 ) 



x |(A^ 2 - A^ 2 ) 2 - -^=(A^~ 2 - A^ 2 )(A^ 3 - A^ 3 ) - (-1)^ 3 (A + - A_) 2 

+ (A^ 2 + A^ 2 ) 

So we obtain 

THEOREM 8 For each j -type QRW, we have 

p$(<p) = ^>l 2 + CjM 2 + mew*?) 

where ip = t [a,(3] G $ and 

1 r 2w 



C A , = ^jJ\apf\e*)+ of \e*)\*d9 

C A , 3 = i- jT (rfV) + crl N \e^))(bA N \e^) + dS™(e»))d0 
1 ^,^n,^.„ _ 1 ^^N )[el e )?de 



1 f27T 



2tt Jo 

with a = b = c = —d = l/y/2, here pj^\z) and r[ N \z) satisfy Eqs. (4-12) and (4-13), and 
C z and E z satisfy Eqs. (4-H) and (4-15). 

Here we consider U = H (the Hadamard walk), (p = and k — 1. From Theorem 8, 

noting that p^\z) = z for any N > 2, we have 

COROLLARY 9 For each j -type QRW, if N > 2, then 

<rV) = K 1 + 27 1" ^"V)! 2 ^) (1 + 2»(^)) 

i N) (e ie )\ 2 d6 





/ 1 r 2 ^ 


\a\ 2 + 


V27T io 



where rf\z) = 0, rf\z) = z 3 / (2 — z 2 ) and for N > 4 



~(AT)/_ A ^ 2 </Af-3(^)^Af-4(^) 



V^-W*)) 2 - -W*) " V^-l)^ 3 
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with 




In particular, when (p — t [0, 1] = \R), k — 1 and N — 2, ... ,6 cases, the above corollary for 
the A-type QRW implies that the values P^i(^[0, 1]) (N = 2, ... ,6) satisfy the conjecture 
given by Eq. (4.6). 
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